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Abstract 


The notion of symmetrization, also known as Davenport’s reflection 
principle, is well known in the area of the discrepancy theory and quasi- 
Monte Carlo (QMC) integration. In this paper we consider applying a 
symmetrization technique to a certain class of QMC point sets called dig¬ 
ital nets over Zb. Although symmetrization has been recognized as a 
geometric technique in the multi-dimensional unit cube, we give another 
look at symmetrization as a geometric technique in a compact totally dis¬ 
connected abelian group with dyadic arithmetic operations. Based on this 
observation we generalize the notion of symmetrization from base 2 to an 
arbitrary base b £ N, b > 2. Subsequently, we study the QMC integration 
error of symmetrized digital nets over Zb in a reproducing kernel Hilbert 
space. The result can be applied to component-by-component construc¬ 
tion or Korobov construction for finding good symmetrized (higher order) 
polynomial lattice rules which achieve high order convergence of the in¬ 
tegration error for smooth integrands at the expense of an exponential 
growth of the number of points with the dimension. Moreover, we con¬ 
sider two-dimensional symmetrized Hammersley point sets in prime base 
b, and prove that the minimum Dick weight is large enough to achieve the 
best possible order of L p discrepancy for all 1 < p < oo. 

Keywords: Quasi-Monte Carlo, 6-adic symmetrization, digital nets, Hammers¬ 
ley point sets 

MSC classifications: 11K38, 65C05. 

1 Introduction 

In this paper we study multivariate integration of functions defined over the 
s-dimensional unit cube. For an integrable function / : [0, l] s -A R, we denote 
the true integral of / by 
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We consider an approximation of 1(f) given in the form 

1 1 X&P 

where P C [0, l] s denotes a finite point set in which we count points according 
to their multiplicity. I(f;P) is called a quasi-Monte Carlo (QMC) integration 
rule of / over P. Obviously the absolute error |/(/; P) — 1(f) | depends only on a 
point set P for a given /. There are two prominent classes of point sets: digital 
nets mm and integration lattices mm- We are concerned with digital nets 
in this paper. 

The quality of a point set has been often measured by the so-called dis¬ 
crepancy mm- For t = (ti,... ,t s ) £ [0, l] s , we denote by [0, t) the anchored 
axis-parallel rectangle [0, t\) x [0, tf) x • ■ ■ x [0, t s ). The local discrepancy function 
A(-; P) : [0, l] s -> R is defined by 

A(t; P) = —— l[o,t)( a: ) — ]^[ tj) 

1 1 xeP j =l 


where l[o,t) denotes the characteristic function of [0, t). Now for 1 < p < oo the 
L p discrepancy of P is defined as the L p -norm of A(-; P), i.e., 


L p (P) := 



i /p 


|A(t; P)| p dt 


with the obvious modification for p = oo. When / has bounded variation V(f) 
on [0, l] s in the sense of Hardy and Krause, the absolute error is bounded by 


|J(/;P)-J(/)|<y(/)i 00 (P). 

This inequality is called the Koksma-Hlawka inequality [TT 1 , Chapter 2], An 
inequality of similar type also holds for 1 < p < oo, see for instance [Uj- This 
is why we consider the L p discrepancy as a quality measure of a point set. 

It is, however, not an easy task to construct point sets with low L p discrep¬ 
ancy. As an example, let us consider the two-dimensional Hammersley point 
sets in base b defined as follows. 


Definition 1. Let b > 2 be a positive integer. For m £ N, the two-dimensional 
Hammersley point set in base b consisting of b m points is defined as 



It is known that Ph has optimal order of the Loo discrepancy, while it does 
not have optimal order of the L p discrepancy for all p £ [ 1, oo), see for instance 
m- Here the general lower bounds on the L p discrepancy for all s £ N have 
been shown by Roth m for p > 2, Schmidt E2H0 for p > 1 and Halasz nni 
for p = 1. Note that the optimal orders of the L\ and discrepancy are still 
unknown for s > 3. 

There are several ways to modify Ph such that the modified point set has 
optimal order of L p discrepancy for p £ [l,oo), see for instance m h m. 
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The notion of symmetrization (also known as Davenport’s reflection principle) 
[2| is one of the best-known remedies in order to achieve optimal order of L p 
discrepancy and has been thoroughly studied in the literature, see for instance 
[nun mm uni- We also refer to 0 in which symmetrized point sets are 
studied in the context of QMC rules using integration lattices. 

Although the original symmetrization due to Davenport has been recog¬ 
nized as a geometric technique in the s-dimensional unit cube, in this paper, we 
give another look at symmetrization as that in a compact totally disconnected 
abelian group with dyadic arithmetic operations. This implies that symmetriza¬ 
tion fits quite well with dyadic structure of point sets and also with the tools 
used for analyzing the point sets. Therefore, it is quite reasonable to consider 
symmetrized two-dimensional Hammersley point sets in base 2 mm , or more 
generally, symmetrized digital nets in base 2 mi- Although there are some 
exceptions as in moss where symmetrization is shown to be helpful even if 
point sets have not dyadic structure, it must be interesting to find a geometric 
symmetrization technique which acts on a compact totally disconnected abelian 
group with 6-dic arithmetic operations for 6 > 2. 

The aim of this paper is two-fold: to generalize the notion of symmetriza¬ 
tion from base 2 to an arbitrary base b £ N, b > 2 and to obtain some basic 
results on the QMC integration error of symmetrized digital nets in base b. In 
particular, we study the worst-case error of symmetrized digital nets in base 
b in a reproducing kernel Hilbert space (RKHS) in Section Q] This result for 
digital nets can be regarded as an analog of the result for lattice rules obtained 
in [SJ Section 4.2] where only the sum of the half-period cosine space and the 
Korobov space is considered as a RKHS. In Section [H we also study the mean 
square worst-case error with respect to a random digital shift in a RKHS. Fur¬ 
thermore, in Section [5j we prove that symmetrized Hammersley point sets in 
base b achieve the best possible order of L p discrepancy for all 1 < p < oo. 

Notation. Let N be the set of positive integers and No := N U {0}. Let 
C be the set of all complex numbers. For a positive integer b > 2, Z;, denotes 
the residue class ring modulo b , which is identified with the set {0,1,..., 6 — 1} 
equipped with addition and multiplication modulo b. For x £ [0,1], its 6-adic 
expansion x = with £ Zb is unique in the sense that infinitely 

many of the Q are different from 6 — 1 except for the endpoint x = 1 for which 
all £i’s are equal to 6—1. Note that for 1 £ N we use the 6-adic expansion 1-6°, 
whereas for 1 £ [0,1] we use the 6-adic expansion (6 — 1)(6 _1 + b~ 2 + ...). It 
will be clear from the context which expansion we use. 


2 Preliminaries 

Here we recall necessary background and notation, including infinite direct prod¬ 
ucts of Zf,, Walsh functions and digital nets (with infinite digit expansions) over 
Z{,. We essentially follow the exposition of [9] Section 2]. 

2.1 Infinite direct products of Z& 

Let us define G = (Zb) N , which is a compact totally disconnected abelian group 
with the product topology where Zf, is considered to be a discrete group. We 
denote by ® and 0 addition and subtraction in G, respectively. Let v be 
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the product measure on G inherited from Z5, that is, for every cylinder set 


E = nr=i %i x rij>n+i 'Lb With Z; C Zf, for 1 < * < n, we have v{E) = 



A character on G is a continuous group homomorphism from G to {z £ C : 
|z| = 1}, which is a multiplicative group of complex numbers whose absolute 
value is 1. We define the fc-th character as follows. 

Definition 2. For a positive integer b > 2, let ui := exp(27r\/—1/6) be the 
primitive b-th root of unity. Let z = (£ 1 , ( 2 , • ■ •) £ G and k £ No whose b-adic 
expansion is given by k = kq + K\b + • • • + K a _i& a_1 with kq, ..., K a -i £ 'Lb- 
Then the k-th character W k : G —> {1 ,u >,... ,w b_1 } is defined as 


W k {z) : = w «°Ci+-+«a-iC 


(1) 


We note that every character on G is equal to some W k , see US- 

Let us now consider the higher-dimensional case. Let G s denote the s-ary 
Cartesian product of G. Note that G s is also a compact totally disconnected 
abelian group with the product topology. The operators © and © now denote 
addition and subtraction in G s , respectively. We denote by v the product 
measure on G s inherited from v. The fc-th character on G s can be defined as 
follows. 

Definition 3. For a positive integer b > 2 and a dimension s £ N, let z = 
(zi,...,z s ) £ G s and k = (ki,...,k s ) £ Ng. Then the k-th character W k : 
G s —> {1, Ub ,..., w£ -1 } is defined as 


S 


W k (z) -.= \{W kj {z 0 ). 


We note again that every character on G s is equal to some W k as with the 
one-dimensional case. 

The group G is related to the unit interval [0,1] through the following maps 
7r : G —> [0,1] and a : [0,1] — > G. Let z = (Ci, C2, • ■ •) £ G, and let x £ 
[0,1] with its unique 6-adic expansion x = £ib~ l with £ Z&. Then the 

projection map 7 r is defined as tt(z) := C ib~ l and the section map a is 

defined as a(x) := (£ 1 ,^ 2 , ■ ■ • By definition, 7r is surjective and er is injective. 
For the s-dimensional case, both the projection and section maps are applied 
componentwise, through which the group G s is related to the unit cube [0, l] s . 
We note that 7r is a continuous map and that n o a = idropp. We summarize 
some important facts below [9j Lemma 4, Propositions 3 and 5]. 

Proposition 4. The following holds true: 

1. For k £ No, we have 



2. For k,l £ Ng, we have 
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3. For f £ L 1 (G S ), we have 



f(z)du(z) 


4- For f £ i 1 ([0, l] s ), we have 



f(x)dx 



f(a(x))dx. 


I f(n(z))dis(z). 

G‘ 


5. Let iJ„ := { 2 : = (Ci, C 2) • ■ ■) € G : Cl = C2 = ■ ■ ■ = ( n = 0}. Then we have 


E 

ke % 
kj<b n ,Vj 


ifzeH s n , 
0 otherwise. 


2.2 Walsh functions 


Walsh functions play a central role in the analysis of digital nets. We refer to 
[0 Appendix A] for general information on Walsh functions in the context of 
numerical integration. We first give the definition for the one-dimensional case. 

Definition 5. For a positive integer b > 2, let 0Jb = exp(27rv / —1/6). We 
denote the b-adic expansion of k £ No by k = kq + k± b + • • • + K Q _i6 a_1 
with Ko, ■ • •, K a -i £ Z;,. TTien f/ie k-th b-adic Walsh function &walfc : [0,1] —> 
{1, u>b, • ■ •, is defined as 

t,wal fc (x) := w” o€i+ '" +k “-i^, 

/or x £ [0,1] with its unique b-adic expansion x = £ 16 _1 + £ 2 b~ 2 + • • •. 

This definition can be generalized to the higher-dimensional case. 

Definition 6. For a positive integer b > 2 and a dimension s £ N, let x = 
(xi,... ,x s ) £ [0, l] s and k = (ki,... ,k s ) £ N§. Then the k-th b-adic Walsh 
function &walfc : [0, l] s —> {l,u>b, ■ ■ ■, u^ -1 } is defined as 

S 

b w&lk(x) := bwalfe^Xj). 

3 =1 

Since we shall always use Walsh functions in a fixed base 6, we omit the 
subscript and simply write walfc or walfc in this paper. From the definitions of 
characters on G s and Walsh functions, we see that for any x £ [0, l] s 

walfc (®) = Wfc(cr(a:)). (2) 

Since the system {walfc : fc £ Nq} is a complete orthonormal system in 
L 2 ([0,1] s ) [2] Theorem A.If], we have a Walsh series expansion for / £ L 2 QO, l] s ) 

E /(fc)walfc, 

where the fc-th Walsh coefficient is given by 

f(k)= [ f{x) walfc(x)dx. 

A[o,i] s 

We refer to El Appendix A.3] and |H1 Lemma 18] for a discussion about the 
pointwise absolute convergence of the Walsh series to /. 
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2.3 Digital nets 

Here we define digital nets in G s by using infinite-column generating matrices, 
i.e., generating matrices whose each column can contain infinitely many non-zero 
entries. This definition has been recently given in [9]. 

Definition 7. For m,s £ N, let C\ ■■ C s £ Z^ xm . For 0 < n < b m , de¬ 
note the b-adic expansion of n by n = Y^iiLo Vib 1 with Vi £ Zf,. Put z n = 
(z n , i, • • ■ ,4, s ) € G s with 

Z'n.j — (^Oi ^1) ■ ■ * i Vm— l) ' Gj , 

for 1 < j < s. Then the set V = {zq, ... , zi,m_ 1 } c G s is called a digital net 
over Z;, in G s with generating matrices C \,..., C s . 

Moreover, the set P := {7r(z): z £ V} C [0, l] s is called a digital net over 
Zh in [0, l] s with generating matrices C \,..., C s . 

We note that every digital net in G s is a Z^-module of G s as well as a 
subgroup of G s . If each column of generating matrices consists of only finitely 
many non-zero entries, the above definition of a digital net over Zf, in [0, l] s 
reduces to that given by Niederreiter El. 

The dual net of a digital net plays an important role in the subsequent 
analysis. For a digital net V in G s , we denote its dual net by V 1 - C Ng which 
is defined as follows. 


Definition 8. Let V be a digital net in G s . The dual net ofV is defined as 

V 1 - := {k = (k u ..., k s ) £ Ng : hC! © • • • ® k s C s = (0,..., 0) e Z™}, 


where we write kj = (ftyo, Kj, i, • ■ •) for kj with its b-adic expansion kj = + 

Kj^b + ■ ■ ■, which is indeed a finite expansion. 

Since Wifs are characters on G s , the following lemma can be established 
from Definition [8l which connects a digital net in G s with characters. 

Lemma 9. Let V be a digital net in G s and V 1 - be its dual net. For k £ Ng, 
we have 


E w ^ z ) 


\V\ ifk£V ± , 

0 otherwise. 


3 The 6-adic symmetrization 

In this section we generalize the notion of symmetrization from base 2 to an 
arbitrary base b £ N, b > 2. Before that, we recall the original symmetrization 
introduced by Davenport [2| and give another look at it as a geometric technique 
in G s with 6 = 2. Let P C [0, l] 2 be a finite point set. Then the symmetrized 
point set in the sense of Davenport is defined as 

psym.D := u (a;, 1 - y) : (x, y) £ P}. 

It is often the case that the symmetrized point set is defined as 

P sym ■= {(a, y) U 0,1 — y) U (1 — x, y) U (1 — x, 1 — y) : ( x , y) £ P}. 
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In the remainder of this paper we only consider the symmetrized point set 
defined in the latter sense. For the higher-dimensional case, the symmetrized 
point set of P C [0, l] s is defined as 

-P sym := (symja:) : x £ P, u C {1,..., s}} , 

where sym u (£c) denotes the s-dimensional vector whose j-th coordinate is 1 — Xj 
if j £ u and Xj otherwise, that is, sym u (a:) = (y i,..., y s ) with 


Vj = 



if j £ u, 
otherwise. 


By definition, we have |P sym | = 2 S |P|. 

Here we give another look at the original symmetrization. Let 6 = 2 and 
z = (Ci, C2, - - -) £ G with Q G Z 2 . By denoting e = (1,1,...) G G, we have 
z ® e = (1 — ( 1 ,1 — ( 2 , ■ ■ •) £ G and thus 


7 t(z © e) 


1 — Ci , 1 — C2 


2 2 2 




C2 

2 2 


1 — 7 r(z). 


For z £ G 5 , we denote by sym(f(z) the s-dimensional vector whose j-th coor¬ 
dinate is Zj © e if j £ u and Zj otherwise. Then the symmetrized point set of 
PC G s can be given by 

P^ m := {sym^(z) :zGP,uC s}}. 

As a natural extension from b = 2 to an arbitrary positive integer b > 2, 
we now introduce the notion of the 6-adic symmetrization. For l £ Zb, we 
write ei = (1,1,...) G G. For a vector l = (li,...,l a ) G Z s h , we write ej = 

( e ii, • • • 1 e h) £ 

Definition 10. For a point set V C G s , its symmetrized point set is defined as 


P sym :={z®ei-.z£V,l£Z s b }. 

For a point set P C [0, l] s , its symmetrized point set is defined as 

P sym = 7r [(<r(P)) sym ] ■ 

By definition, we have |P sym | = b s \V\ and |P sym | = 6 S |P|. 


3.1 Symmetrized digital nets 

In the following let P be a digital net in G s . We write Q = (e* : l G Z b }, 
which is also a digital net in G s . From Definition 1101 the symmetrized point set 
P sym can be regarded as the direct sum of two digital nets P and Q. Thus, it 
is obvious that the following holds true. 

Lemma 11. For a digital net P in G s , let , p syra be the symmetrized point set 
of P. Then P s v m is also a digital net in G s . 
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Remark 12. Consider a digital net V in G s constructed with infinite-row gen¬ 
erating matrices C \,..., C s G Z^ xm . Then 'P sym is a digital net in G s whose 


generating matrices D± 



are 

given 

as 





1 ... j 

- 1 

j j 

+ 1 • 

s 




(0 ••• 

0 

1 

0 

■ o\ 




0 ••• 

0 

1 

0 

• 0 

G Z^ xs 

Dj = (Cj,Ej) with 

Ej = 

0 ••• 

0 

1 

0 

• 0 







• ) 



For a point set 'P sym , we have the following orthogonal property. 


Lemma 13. Let V be a digital net in G s and , P sym be its symmetrized point 
set. For k G Ng, we have 


E w ^ z ) 

Z^.'psym 


|P sym | ifkGV ± (l£ s , 
0 otherwise. 


In the above, £ := {k G No : S(k) = 0 (mod b)} where 6{k) denotes the b-adic 
sum-of-digits of k and is given as S(k) := kq + Ki + ■ • • for k = n o + n\b + • • • . 

Proof. From the definition of the 6-adic symmetrization, we have 


E ^( z ) = EE fffe ( z0e! ) 

zg -psym ZGVIGZI 

= E Wk(z ^ E 

zG V l&l 


Since V is a digital net in G s , the first sum in the last expression equals \V\ if 
k G V and 0 otherwise. On the second sum in the last expression, we have 

E Wik(ei) = EII4 i%) = n E u b S{ki) - 

j —1 1 Ij 

As cob denotes the primitive 6-th root of unity, the inner sum on the rightmost- 
side above equals b if 6(kj) = 0 (mod b) and 0 otherwise. Thus, 


E Wfc(ei) 

All together, we obtain 


b s if k G £ s , 
0 otherwise. 


E w ^ z ) 

z£V sym 


b s \p\ if keP ± n£ s , 

0 otherwise. 


Since |P sym | = b s \P\, the result follows. 


Combining Lemmas [9] [TT] and [13] implies that 


(P sym )-i- = p ± n£ s . 
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Remark 14. In fact, the argument in this subsection can be generalized in the 
following way. Let V and V be digital nets in G s . Consider the direct sum 

1Z= {z® z' : z £V,z' £ V'}. 


Then IZ is also a digital net in G s and satisfies the orthogonal property 


E w ^ z ) 

z&TZ 


\n\ ifkGV ± n(V) ± , 

0 otherwise. 


In the remainder of this paper, however, we only consider the case V' 
which gives us TZ = p s y m . 


Q, 


4 QMC integration over symmetrized digital nets 


Let us consider a reproducing kernel Hilbert space (RKHS) H with reproducing 
kernel K : [0, l] s x [0, l] s —y R. The inner product in % is denoted by ( f,g)u for 
f,g €% and the associated norm is denoted by WfWn '■= \/ (/, f)u- It is known 
that if Jj Q -n, y/K(x, x) dec < oo the squared worst-case error in the space H of 
a QMC integration over a point set P C [0, l] s is given by 


e 2 (P,K) := 


sup |/(/)-/(/;P)| 
\ll/llw<l 


I K( Xl y)dxdy- E / K{x,y) dy + tE E K ( x >y) 
J l ° d] 2a 1^1“ x.wgp 


(3) 


see for instance [25]. For k,l £ Ng, the (k, Z)-th Walsh coefficient of K is defined 
by 


K(k,l) 



K(x, y)walfc(a;)walj (y) da; dy. 


In the following we always assume Jj Q 1 , 3 y/K(x, x) da; < oo. We study the 
worst-case error of symmetrized point sets in a RKHS first and then move on 
to the mean square worst-case error with respect to a random digital shift. 


4.1 The worst-case error 

From the proof of [jJJ Proposition 21], we have the following pointwise absolute 
convergence of the Walsh series of K. 

Lemma 15. Let K be a continuous reproducing kernel. We assume 

E i£(m)i<«>. 

k,ie% 

For any z,w G G s , we have 

K(tt(z),tt(w)) = E K{k,l)W h (z)Wi(w). (4) 

fc,2e Ng 
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Under some assumptions on K, the worst-case error is given as follows. 

Theorem 16. Let 'P,'P _L be a digital net in G s and its dual net, respectively, 
and V sym be the symmetrized point set of V . Let K be a continuous reproducing 
kernel which satisfies 

/ \/ K(x,x)dx < oo and \K(k,l)\<oo. 

MeNg 

The squared worst-case error of a QMC integration over 7r('P sym ) is given by 

e 2 (Tr{V syin ),K) = J2 k (M)- 

fe,ie7’ J -n£»\{o} 

Although the proof is almost the same with that of !9j Proposition 21], we 
provide it below for the sake of completeness. 

Proof. We evaluate the three terms of ([3]) separately in which n(P sym ) is sub¬ 
stituted into P. The first term of ([3]) is simply 


'[o.il 2 


K(x,y) dx d y = K( 0,0), 


by the definition of the Walsh coefficients. For the second term of ([3]), we have 
2 


\V S 


1 


\P 


sym 


\v s 


E 

zG'psym 

E 

£-psyi 

■ E 

Z^.'psym • 


'[0,1] S 


K(ir(z),y)dy 


1 


zcV s ym 


W K[ ” (z) ' v)Av+ V-y- 

/ K(tt(z), tt(w)) dv(w) + 


E 

ZgP" jm 

1 


IP 


sym 


•Wl" 

E 


K{y,n(z))dy 


K(tt(w),tt(z)) du(w) 


z£V s ym 1 


I V 


sym 


E E K(k, l)Wk(z) / Wi{w) du{w) 


| 7 ?sym 

= y 


z6T>'f” k,lG Ng ^ 

E E K(k,l)W t {z) f Wfc(tc) du(w) 

z £V syrn k,lG% JG‘ 

1 


fceNg 


\V S 


E w.W + EUo.Djpss, E WM 


z£V sym 


iem 


z£V sym 


= Y K(k,o)+ y 

fee'P-'-nf 8 «e'P- L n£ s 

where we use the symmetry of K, Item 4 of Proposition [4] Equation (|4]), Item 1 
of Proposition []] and Lemma [13] in this order for each of the five equalities. 
Finally, for the last term of f3]), we have 

ypi ^12 E K(TT(z),n(w)) 

1 1 2,«>G? sym 

= ip^p E E K(k,i)w*(zjmiw) 

1 1 2,uiE'P sym k,le% 
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= £ 

k,ie Ng 1 1 zev sym 1 1 wev sym 

= E Em), 

fc,ieP J -nf s 

where we use Equation Q and Lemma m in the first and third equalities, 
respectively. Substituting these results into the right-hand side of © , the result 
follows. □ 

Remark 17. The result of Theorem \16\ has some similarity to that o/© Propo¬ 
sition 22]. When a QMC integration over a folded digital net by means of the 
b-adic tent transformation is considered, the squared worst-case error in a RKHS 
is given by 


E El k/b\,n/b\), 

fe,JG'P- L nf s \{o} 

where we write |ycj = (LaqJ,. .., |_£ S J) for x = {x\,...,x s ) € R s . Since the 
number of points of a folded digital net is the same as that of an original digital 
net, a folded digital net has a cost advantage over a symmetrized digital net. 

4.2 The mean square worst-case error 

Here we study the mean square worst-case error of symmetrized point sets with 
respect to a random digital shift. Now the error criterion is given by 

e 2 ( 7T (P sym ),K)= [ e 2 (7r('P sym ) © er, K) der, 

•'[0,11 s 

where the operator © is defined for x, y £ [0, l] 5 as 

x © y := 7T [a(x) © a (y)\. 

From [6] Theorem 12.4], we have the following. 

Lemma 18. For a point set V C G s and a reproducing kernel I\ £ L2([0, l] 2s ), 
the mean square worst-case error of a set tt(V) with respect to a random digital 
shift is given by 


e 2 ^(V),K)=e 2 ^(V),K ds ), 


where K ds is called a digital shift invariant kernel defined as 


K ds (x,y) 



K(x © cr, y © cr) der, 


for any x, y £ [0, l] s . 

Similarly to Lemma fTSl we have the following pointwise absolute convergence 
of the Walsh series of K ds . 
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Lemma 19. Let K be a continuous reproducing kernel. We assume 

y \K(k, k )| < oo. 

fcSNg 


For any z,w £ G s , we have 

K ds (n(z),n{w)) = ^2 K(k,k)W k (z)W k (w). (5) 

feGNg 

Proof. By the assumption |^(fe)fe)| < oo, the right-hand side on © 

converges absolutely. Thus it suffices to show that 

lim Y' K(k,k)W k (z)W k (w) = K da (-K(z),n(w)). 

n—t oo z ' 
ke% 
kj <b n ,Vj 


In fact we have 


Y2 K(k,k)W k {z)W k (w) 
ke Nq 

kj <b n ,vj 

= Y2 W k {z)w k (w) 

ke No 
kj<b n ,Vj 

= y w k ( Z )W^) 

fcG Ng 
kj<b n ,Vj 

K(n(z'),- tt(w')) Y2 W k {{z Q z') Q (w Q w'))du{z')Au{w') 
fee Nq 

kj<b n ,Vj 



/G 2s 


K(tt(z '), 7r(m , ))W / fc(;s')IPfc(u/) di'(z') dt'(tu') 


/[ 0 , 1] 2 


Iv (a;, y)wal fc (a;)wal fe (y) da; dy 


/G 2s 


K(tt(z © z r ), ir(w © w')) 




W k (w' © z') d u(z') &u(w'), 


kj<b n ,Vj 


where we use Item 4 of Proposition [4] twice in the second equality. Similarly to 
the proof of [5J Proposition 21], let us define two sets H n := {z = (Ci 5 C 2 , • ■ ■) £ 
G : Ci = C 2 = • • • = (n = 0} and J„ j2s := {( z , w) £ G 2s : w © z £ H^}. Then 
from Item 5 of Proposition |4] we have 


Y2 W k (w'Qz') 

kj <b n ,Vj 


b sn if (z',w')£ J n , 2 s, 
0 otherwise. 


Thus we have 


Y2 K(k,k)W k {z)W k (w) 
ke Nq 

kj <b n ,\/j 

= b sn f K(n(z © z'), 7 t(w © w')) du(z') du(w') 

J Jn,2s 
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K(tt(z ® z'), n(w ® «/)) d v(z') du(w') 


f 6» /' 

JG a Jw'QH a 


—>■ 



K(n(z ® w'), 7 t(w ® w')) dv(w') 


as n —> oo, 


where the last convergence stems from the facts that K o 7r is continuous since 
both K and 7r are continuous and that the product measure of the set w' QH° n 
equals b sn for any w' £ G s . Finally we have 


lim y K(k,k)W k (z)W k (w) 

n—too z —* 

ke% 

kj<b n ,\/j 


= / K(n(z ® w'), n(w ® w')) du(w’) 

Jg s 

= / K(tt(z) ® ir(w'), tt(w) ® n(w')) du(w l ) 

Jg » 

= / K(tt(z) ® ® cr)dcr = K ds (Tr(z) 1 TT(w)), 


where we use Item 3 of Proposition U in the third equality. Thus the result 
follows. □ 


Under some assumptions on IF, the mean square worst-case error with re¬ 
spect to a random digital shift is given as follows. 

Theorem 20. Let V,V _L be a digital net in G s and its dual net, respectively, 
and 'p sym be the symmetrized point set ofV. Let K be a continuous reproducing 
kernel which satisfies 

\JK(x , x) dx < oo and | K(k, k) | < oo. 

fee Ng 

The mean square worst-case error of a QMC integration over 7r('P sym ) with 
respect to a random digital shift is given by 

e 2 (ir(V sym ),K)= K(k,k). 

fceP- L n£' , \{o} 

Proof. We evaluate the three terms of m separately in which 7r('P sym ) and K ds 
are substituted into P and K, respectively. The first term of Q is given by 




K ds (x,y) da: dy 


' G 2a 


K ds (ir(z), 7r(w)) di j(z) du(w) 


Y, K(k,k ) [ W k (z)W k {w)di>(z)di/(w) 
J G 2s 


fee ns 


= IF(0, 0), 

where we use Item 1 of Proposition 0] in the third equality. For the second term 
of dH), we have 

E [ u K Mz),y)dy=-^- y f K ds (7r(z),7r(w))diy(w) 
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= |-p sym | E E K{k,k)W k (z) w k (w) du(w) 

I I zg -psym fcgpjg ^ G “ 

= 2K(0, 0), 

where we use Item 4 of Proposition © Equation © and Item 1 of Proposition 0] 
in this order for each of three equalities. Finally, for the last term of Q, we 
have 

Tps Li E K ds (n(z),n(w)) 

1 1 z,wev sym 

= yp^V 2 E E K(k,k)W k (z)W k (w) 

1 1 z.-iue-p^ 111 fee Ng 

= K(k,k )—53 ^(z)—53 W k (w) 

ke% ' ' zev sym 1 1 wev sym 

= 53 K(k,k), 

kev ± n£ B 

where we use Equation © and Lemma m in the first and third equalities, 
respectively. Substituting these results into the right-hand side of © , the result 
follows. □ 

Remark 21. The result of Theorem \20\ has some similarity to that of © The¬ 
orem 16]. When a QMC integration over a “digitally shifted and then folded” 
digital net is considered, the mean square worst-case error with respect to a 
random digital shift in a RKHS is given by 

53 K([k/b\, \_k/b\), 

feeP- L nf s \{o} 

where again we write \x\ = ([xij,..., |yc«J) for x = (xi,...,x s ) £ R s . As 
already mentioned in Remark m this approach has a cost advantage over a 
symmetrized digital net. 

Remark 22. Following the same arguments as in © Sections 4-6], Theo- 
rems \lb\ and \20\ can be applied to component-by-component (CBC) construction 
or Korobov construction of good symmetrized (higher order) polynomial lattice 
rules which achieve high order convergence of the worst-case error in an unan¬ 
chored Sobolev space of smoothness a £ N, a > 2. For instance, for m £ N and 
a prime b, the CBC construction can find symmetrized higher order polynomial 
lattice rules with b m+s points which achieve the worst-case error convergence of 
0(b~ arn+e ) with an arbitrary small e > 0. Moreover, this construction can be 
done in 0(samb am / 2 ) arithmetic operations using 0(b am / 2 ) memory. 

5 Discrepancy bounds for symmetrized Ham- 
mersley point sets 

Finally in this paper we prove that symmetrized Hammersley point sets in 
prime base b achieve the best possible order of L p discrepancy for all 1 < 
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p < oo. According to Definitions [T] and 1101 the two-dimensional symmetrized 
Hammersley point set in base b is given as follows. 

Definition 23. Let b >2 be a positive integer. Form £ N, the two-dimensional 
symmetrized Hammersley point set in base b is a point set consisting of b m+2 
defined as 

P H m '■= {( x o.,y a ) ■■ a = (a i,... ,a m+2 ) £ {0,1,.. .,6- l} m + 2 } , 


where 


! _ oi®a m _|_i . . a m ®a m _|_i . a m _|_i . a m _)_i , 

J a — 5 \ * * * I b m I frm+1 I b m + 2 ' ' ' ‘ 5 

_ a m ®a m _)_2 I I «l®a m -|-2 I «m + 2 I CLm -(-2 

y<x — - 5 -r • • • i -r 6 m+i -I- T-• 

We also consider a truncated version of Pff m . For a positive integer n > 
m + 2, the two-dimensional truncated symmetrized Hammersley point set in 
base b , denoted by Pff m,n , is given as 

P™ := {«, v2) : a = (a ls ..., a m+2 ) £ {0,1,..., b - l} m+2 } , 


where 


„n _ fll8flm+l | ... I a m®Qm+l _i_ Qm +1 i i Q.m+1 

•' a 5 r * ' * i - frm i j,m +1 r * ' * “r , 

n Qm®Qm + 2 | I , Qm + 2 ■ ■ «m+2 

Ua b r * * ' i 5 m ' 6 m + 1 i * * ‘ i 5 n 


Here note that Pj^ m ’" is a digital net over Zf, with generating matrices of size 
n x (m + 2) 


/1 

0 

... 0 

1 

0 ^ 


( 0 

... 0 

1 

0 

1 \ 

0 

1 

... 0 

1 

0 


0 

... 1 

0 

0 

1 

0 

0 

... 1 

1 

0 

;C 2 = 

1 

... 0 

0 

0 

1 

0 

0 

... 0 

1 

0 


0 

... 0 

0 

0 

1 

V o 

0 

... 0 

1 

0 ) 


V o 

... 0 

0 

0 

1 ) 


( 6 ) 


The following theorem is the main result of this section. 


Theorem 24. For a prime b and m £ N, let Pff m be the symmetrized Ham¬ 
mersley point set in base b consisting of N = b m+2 points. Then, for all 
1 < p < oo, the L p discrepancy of P^ m is of the best possible order. Namely, 
for any 1 < p < oo there exists a constant C p > 0 such that 


L P (P s H n < c p 


yj TYl + 2 
bm +2 


si V lo & b N 

Cp N 


Our proof consists of two parts. We first prove that the L p discrepancy of 
psym.n - g possible order for any 1 < p < oo when n > 2m. Then 

we show that the difference between the L p discrepancies of Pff m and pff m,n is 
small enough that the L p discrepancy of Pff m is still of the best possible order 
for any 1 < p < oo. 
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5.1 The minimum Dick weight of point sets 

To prove the first part, we introduce the Dick weight function /i 2 given in 001 
and the minimum Dick weight p 2 (P) for a two-dimensional digital net P. 

Definition 25. For k € N, we denote its b-adic expansion by k = n\ 6 ai_1 + 
K 2 b a2 ~ l + • • • + n v b av ~ 1 such that ki, ... ,k v € {1, — 1} and ai > 02 > 

• • • > a v > 0. Then the Dick weight function ^ 2 : No —>■ R is defined as 

! ai + a 2 if v > 2, 

ai ifv = 1, 

0 if k = 0. 

For vectors k = (fci,/c 2 ) £ Nq, we define P2 (fe) := ^ 2 (fci) + ^2(^2)- Moreover, 
let P be a two-dimensional digital net over Zb- Then the minimum Dick weight 
P 2 {P) is defined as 


P 2 {P) ■= min /i 2 (fc)- 
fcepJ-Uf 0,0)} 

The following lemma shows how the minimum Dick weight of a digital net 
connects with a structure of its generating matrices, see [5] Chapter 15]. 

Lemma 26. For m, n £ N with n > m, let P be a digital net over Zb with 
generating matrices C\, C 2 of size n x to. For j = 1,2 and 1 < l < n, let c^i 
denote the l-th row vector ofCj. When l > n, Cy; denotes the vector consisting 
of m zeros. Let p be a positive integer such that for all 1 < ii tVl < ■ ■ ■ < h, 1 < 
2to and 1 < i 2 ,v 2 < ■ ■ ■ < * 2,1 < 2 m with Vj G No and 

min(t;i, 2 ) min(t>2,2) 

- P ’ 

;=i z=i 

t/ie vectors ci i j l vi ,..., c 2) i 2> „ 2 ,..., c 2 ,i 21 are linearly independent over 

Z;,. T/ien we ftoue p 2 (T > ) > p. 

Although Definition [25] and Lemma [26] are considered only for the two- 
dimensional case, they can be generalized to the s-dimensional case with any 
s G N. Recently, Dick [4] proved that digital nets over Z 2 with large minimum 
Dick weight achieve the best possible order of the L p discrepancy for 1 < p < 00 
and for any number of dimensions. For the two-dimensional case, his result 
also implies the best possible order of the L\ discrepancy with respect to the 
general lower bound by Halasz m ■ Dick’s result was generalized more recently 
in [T5j to digital nets over Zb for a prime b. We specialize their results m 
Corollary 2.2] and [151 Corollary 1.8]) on the L p discrepancy of digital nets for 
the two-dimensional case. 

Proposition 27. Let P be a digital net over Zb consisting of b m points which 
satisfies p 2 (P) > 2 m — t for some integer 0 < t < 2m. Then for all 1 < p < 00 
there exists a constant C p ,t which depends only on t and p such that we have 


m 


L p {P) < C Ptt b m ■ 
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From Lemma [M] and Proposition [23 as the first part of the proof of Theo¬ 
rem [Ml it suffices to show the linear independence properties of the generating 
matrices © such that the minimum Dick weight of pff m,n is large. 

Lemma 28. For m,n £ N with n > 2m, let C\ and C 2 denote the generating 
matrices W- For j = 1,2 and 1 < l < n, Cj t i denotes the l-th row vector of Cj. 
The following sets of the vectors are linearly independent over Z;,: 


1 . 


' , Ci, r> C2,l, • • 

. ,c 2jm+ i_ r } for 0 < r 

< m + 

1, 


2. 


5 Cl,m+l ? C2,r 

} CL^d {C2,l ? • • • 9 C 2 ,rn+l 9 Cl,7*} 

for 1 < r 

< 777 , 

3. 


, , Ci >r , C2,l, • • 

• ) C 2 ^m—n Cj.} for j — 

1,2, 0 

< r < m, 

and ?77 + 1 < 


s < n, 






4 - 


9 c l,ri,2 5 

1, C 2 ,i, . .., c 2>r . 2 2 ,c 2jr . 2 

1} for 

0 < r 1)2 

VI 

s~ 

V 


and 0 < 

^2,2 < r 2 ,i < 

to such that rip + ri,; 

2 + r 2 ,1 

+ r 2 ,2 < 

2 m + 1. 


Since the proof is quite similar to that of [3 Lemma 3.1], we omit it. Using the 
linear independence properties of © shown in the above lemma together with 
Lemma [M] and Proposition [23 we have the following. 

Proposition 29. For m,n £ N with n > 2m, let Pff m ’ n be the truncated 
symmetrized Hammersley point set in base b consisting of N = b m+2 points. 
The minimum Dick weight of Pff m ’ n is larger than 2m+ 1, that is, 

P2{PTn >2 m+l. 

This implies that the L p discrepancy of p s ^ m,n is bounded by 

t r psym,n\ ^ \/m + 2 ydogj, N 

^p\^h J 5s o P bm+2 ~^p N = 

for any p £ [l,oo) with a positive constant C p depending only on p. 

Since the proof on the minimum Dick weight of p^ m ’” is quite similar to that of 
[3 Lemma 3.3], we omit it and just give some comments on the L p discrepancy 
bound instead. As pff m ’ n consists of b m+2 points and the minimum Dick weight 
of pff m ’ n is larger than 2m + 1 = 2(m + 2) — 3, the t -value in Proposition [271 of 
jjsyiri.n a i wa y S e q Ua j s 3 f or an y m g pj. This is why we simply write C p instead of 
C Pit in the above proposition. This t-value is same as that for two-dimensional 
folded Hammersley point sets as given in [3 Lemma 3.3]. Moreover, if one wants 
to get an explicit value of C p , it might be better to use the Littlewood-Paley 
inequality in conjunction with the Haar coefficients of the local discrepancy 
function for pff m ’ n a s done in m, instead of our proof based on the linear 
independence properties. However, this is beyond the scope of this paper. 

5.2 Effect of the truncation on the L p discrepancy 

As the second part of the proof of Theorem [Ml we show that the difference 
between the L p discrepancies of and pff m ' n is small enough that the L p 
discrepancy of Pj^ m is still of the best possible order for any 1 < p < oo. Namely 
we show the following. 
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Proposition 30. For m,n £ N with n > m + 2, let P^J m be the symmetrized 
Hammersley point set in base b consisting of N = b m+2 points and p^ m,n be 
its truncated point set. Then we have 

M p r) < WOT"") + (?) 

for any p £ [ 1 , oo). 

Before providing the proof of Proposition [jJOj it should be mentioned that we 
arrive at the result of Theorem [2H by combining Propositions [2H and [HI since 
the second term on the right-hand side of m is small enough for any p £ [ 1 , oo) 
that it does not affect the order of the L p discrepancy when n > 2m. 


Proof. From the definition of the L p discrepancy we have 

l p (pd= (f m-,p s r)\ p dt) 

\J[o, d 2 / 

= ( f | A(i; Pf m ) - A(t; Pj+ A(t; Pj m ’ n )\ p dt) 

Vpm 2 J 

,(PZT' n ) +( f l A (*; p H m ) - A (*; Pf m ’ n )r dt) , 


\ 1/p 


A Ppv 


where we use the Minkowski inequality. Thus, we shall focus on the second term 
on the rightmost-hand side in the following. 

Now for any t £ [0, l ] 2 we have 

A(t; Pff m ) - A (t;Pff m ’ n )=y^ (l[o,t) (x a , 2/a) - 1 (0 , t)«, 2/a)) ■ 

ae{0,l,...,6-l} m + 2 


For the summand on the right-hand side, we have 


l[0,t) (Xai Da) 




l 

-l 

o 


if (x a ,y a ) £ [0,t) and (x^,y^) [0,t), 

if (x a ,y a ) i [0,t) and «,</£) £ [0,t), 
otherwise. 


For a given set a = {ai,..., a m + 2 }, we have 


x a <x a = 


a, 1 © dm+l 


a m ® 


= xl + 


&m+l 
6”(6 — 1 ) 


<x n n + 


Q"m -\-1 ^m+1 &m+l 

^m+1 ^n+1 


Similarly we have 


2/a < 2/a 


&ra +2 

1 ) 



This implies that we never have the case where (x a ,y a ) £ [0,t) and (x™,y™) ^ 
[0,t). Thus, for any t £ [0, l ] 2 we have 


|A(t;Pf m )-A(t;P™)| 
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_ | {a € { 0 , 1, ■ ■ ■, 6 - l} m + 2 : (x a ,y a ) j [0, t) and (x", yg) € [0, t)} \ 

fom+2 

_ I {g € { 0 ,1, ■ ■ ■, b - l} m+2 : x” < ti < x a and y" < t 2 < y a ) \ 

tyn+2 

^ | {a £ { 0 , 1 ,. .., 6 - l } m + 2 : x n a < t x < x” + 1 / 6 ” and y n a <t 2 < y% + 1 / 6 ”} | 
— b m +2 

From the definition of P^ m ’ n , the first coordinate x” takes the value of the form 

ci c 2 {_L_ 1 _ 

b m + b - 1 \b m 6 ” 

for Ci,C 2 £ No with 0 < Ci < b m and 0 < C 2 < 6 . Since x£ does not depend 
on a m + 2 , each value repeats 6 times. The length of the interval between two 
distinct consecutive elements of x™ equals either 

_L('_L-±\ or ±. 

6-1 \b m 6 ") 6 " 

Note that the former length is greater than 1/6” when n > m+2. This also holds 
true for the second coordinate y™. Therefore, when there exist c\,c 2 ,di,d 2 £ No 
with 0 < ci, d\ <b m and 0 < c 2 ,d 2 < b such that 

Cl C 2 ( l 1 \ t < Qt C 2 / 1 1 \ 1 

6 ™ 6—1 V 6 m b n ) 1 “ b m 6 - 1 \b m 6 ") 6 ” 

and 

2 l + d2 (2 _ 2 \ < t < — + d2 (2 _ 2 \ + 2 

b m b-1 \b m 6 ") 2 ~ b m b - 1 \ b m 6 ”) 6 ” ’ 

there are at most 6 2 points of which satisfy x” < ti < x” + 1 / 6 ” and 

Va<t 2 <ya + 1/6". Note that the Lebesgue measure of the set consisting of 
such t is given by 1 /fr 2 ( n_Tn_1 ). On the other hand, when there do not exist 
ci,C 2 ,di,d 2 £ No with 0 < Ci,di < b m and 0 < c 2 ,d 2 < b which satisfy the 
above condition, there is no point of P// rn ’” which satisfy x” < t\ < x” + 1/6” 
and y” < t 2 < y” + 1/6". Therefore, we have 

(L la(t; ^’> ■- a<t; ^ ”'“ )|p dt )^ 1,p 

1 

^ra+ 2 (n—m—l)/p ’ 

Hence the result follows. □ 
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